Empirical Fit to electron-nucleus scattering. 
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An empirical fit to electron-nucleus scattering for A > 2 is made based on world data. It is valid 
for < W < 3.2 GeV and 0.2 < Q 2 < 5 GeV 2 , and can be used with caution at lower Q 2 . The fit 
is based on previous empirical fits to electron-proton and electron-neutron scattering, taking into 
account the effects of Fermi motion plus a substantial extra contribution that fills in the dip between 
the quasi-elastic peak and the A(1232) resonance. 
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There are many applications in nuclear/particle physics for a reliable parametrization of inclusive electron-nucleus 
scattering. One example is the reliable evaluation of radiative corrections to measured data to extract inclusive 
electron scattering cross sections. Another is the evaluation of background in electron scattering from polarized 
ammonia. In this case, the contributions from the unpolarized nitrogen relative to the polarized proton must be taken 
into account to extract spin-dependent inclusive cross section asymmetries. 

A prime motivation for a new empirical fit is the availability of a large body of new, high precision electron-nucleus 
scattering data from Jefferson Lab [lj ■ The present fit was used in the evaluation of the radiative corrections for much 
of these data. 

Our basic fit form is similar to that of a previous fit to inclusive electron-deuteron scattering [3j • 

II. DEFINITIONS AND KINEMATICS 

In terms of the incident electron energy, E, the scattered electron energy, E , and the scattering angle, 9, the 
absolute value of the exchanged 4-momentum squared in electron-nucleon scattering is given by 



Q 2 = (-q) 2 =4EE sin 2 -, 



(1) 



and the mass of the undetected hadronic system is 



% 



W 2 = M 2 



2M p v 



Q 2 



(2) 



with M p the proton mass, v = E — E' , and the small terms involving the electron mass squared have been neglected. 
In the one-photon exchange approximation, the spin-independent cross section for inclusive electron-nucleon scat- 
tering can be expressed in terms of the photon helicity coupling as 



da 



a 2 F cos 2 {9/2) 

dndE' ~ [2Esm 2 {9/2)] 2 



(W 2 (W 2 , Q 2 ) + 2 tan 2 (9/2)W 1 (W 2 , Q 2 )} 



(3) 



where chf — 1/137 is the fine structure constant, and the structure functions W\ and Wi are each the sum of a 
quasi-elastic (W® E ) and an inelastic (W 1 ) piece, as detailed below. The quasi-elastic fit is given in terms of the F\ 
and F 2 structure functions, which are related to W\ and W 2 by F\ = MW\ and F 2 = vW 2 . 
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III. INELASTIC FIT FUNCTION 

The functional form of the fit is: 

W( = (W[(W 2 , Q 2 ) + W? IEC (W 2 ,Q 2 ))f EMC (x') 

where Wf is a Fermi-smeared sum of free proton and neutron contributions, W^ is an additional term to fill in 

the dip between quasi-elastic and A(1232) peaks, and Jemc is an unpublished parametrization of the "EMC" effect 
(nuclear dependence of Deep Inelastic Scattering) by S. Rock in 1994. 
To obtain W 2 , we use the relation: 

Wi = Wl[l + R A (W, Q 2 )}/(1 + v 2 /Q 2 ) 

where 

Ra(W, Q 2 ) = Rp(W, Q 2 )(l + P 6 + P 23 A) 

where R p is the ratio of longitudinal to transverse cross sections from the fit to world proton data given in Ref. [I| ■ 
The fit parameters Pi are listed in Table UlTl 

A. Inelastic term W[{W 2 ,Q 2 ) 

To obtain the inelastic term, Wf(W 2 ,Q 2 ), we use the equation: 

Wf(W 2 ,Q 2 ) = C(x)Y J \ZWl((Wlf 1 Q 2 ) + (A- Z)W^{{W^) 2 1 Q 2 )]f l (4) 

i 

where the correction term depending on x Bjorken is given by: 

C{x) = 1 + P 13 x + P 14 x 2 + P 15 x 3 + P 16 x 4 + P 17 x 5 , (5) 

and Wf and W™ are the free proton [J] and neutron [3| structure functions. The shifted values W[ are given by 

{WD 2 = W 2 + £ik F \q\ - 1E s {y + M) (6) 

where |g| 2 = Q 2 + v 2 , 

U = -3 + 6(t - l)/98 (7) 

and 

fi = 0.0245e ( ~^ /2) . (8) 

The sum is nothing more than a step- wise integration over a Gaussian whose width is controlled by a Fermi momentum 
kp, truncated at ±3er, with a shift in central W related to the binding energy E s . The values of kp and E s used for 
the different nuclei are given in Table Q] 

A k F (GeV) E a (GeV) 

3 0.115 0.001 

3 < A < 8 0.190 0.017 

7 < A < 17 0.228 0.0165 

16 < A < 26 0.230 0.023 

25 < A < 39 0.236 0.018 

38 < A < 56 0.241 0.028 

55 < A < 61 0.241 0.023 

A > 60 0.245 0.018 

TABLE I: Values of Fermi-broadening parameter k F and binding energy parameter E s for different nuclei. 



B. "MEC" term 

To fill in the dip between the quasi-elastic and A(1232) resonance peaks, we added an extra term, which we dubbed 
the "MEC" (meson-exchange current) term. The importance of this term grows with A. The form of W^ is: 

w mec _ ^0_ e -i(VW I ~P 1 ) 2 ]/P2] (g) 

where 

/ = (1 + (Q 2 )')/P 3 f* v p * (1 + P 1S A^ 1+P ^) (10) 

where (Q 2 )') is the larger of 0.3 GeV 2 or Q 2 and x = Q 2 /2Mv. 

C. Jemc 



The function Jemc is given by: 



where 



I EMC = cA a 



a = ao 



E^M* 



where x' is the smaller of 0.7 and Q 2 jIMv and the coefficients on are given in Table HH and where 

„0. 0169+0. 01809a:'+0.0050427(x') 2 



c = e 



The Jemc fit is illustrated for five representative nuclei in Fig. [TJ Note that the original fit was only for data with 
Bjorkcn x < 0.7, that is why we "freeze" the results above that value 



ao &1 Q2 «3 »4 &5 OLfi «7 Ol$ 

-0.069887 2.1888 -24.667 145.29 -497.23 1013.1 -1208.3 775.76 -205.87 
TABLE II: Values of the a t EMC fit parameters. 



IV. FIT PARAMETERS 

The fit parameters are given in Table Hill except for Pig. This parameter, which is important in the strength of the 
MEC term, was fit individually for 3 H e, He, C, Al, and Cu. We assume nearby nuclei will have the same parameter 
values. The results are given in Table HVl 

V. QUASI-ELASTIC CONTRIBUTION 

The quasi-elastic contribution is calculated using the equations in Ref. 5]. The free nucleon factors are taken 
from Ref. [6|]. These form factors are based on inclusive electron scattering, and thus have the 2-photon corrections 
appropriate to the present quasi-elastic fit. The values of Fermi motion parameter kp and binding energy E s are 
slightly different than in Ref. Q: we used those in Table Q] We used a multiplicative Pauli suppression factor given by 

(3/4)(|«1/* F )(l - (|$1/fc,) a )/12) 

for \q\ < 2fci?, otherwise no correction was made. We assumed the same suppression for W\ and Wz- For the scaling 
function, we used Q 

F(V>') = 1.5576/M1 + 1 -7720 2 0// + 0.3014) 2 )(1 + e - 2 - 4291 ^) (11) 
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FIG. 1: Illustration of the ^-dependence of the function Jemc for five values of atomic number A. 



Po 


0.005138 


Pi 
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P, 
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Pie 


-1.905500 


Pl7 
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Pis 


see Table 2 P 19 


-0.045536 


P20 


0.249020 


P21 


-0.137280 


P22 


29.201000 


P 23 


0.004928 



TABLE III: Values of the fit parameters. 



We used the same function for both W\ and Wi. 

The nominal results from the above were corrected as follows: 



Fi = pno m{1 + pj + p sy + p 9y 2 + p wy 3 + p iiyi) 



and 



R = R nom (l + P l2 ) 
where R nom is the ratio of longitudinal to transverse cross sections, defined using: 

pnora = ( j P™ om /l/)(M/F 1 " 0m )(l. + V 2 / ' Q 2 ) - 1 

and where y = (W 2 — M 2 )/\q\. The fit parameters Pj to P\i are given in Table ITTT1 
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TABLE IV: Values of Pis for different nuclei. 



VI. COULOMB CORRECTIONS 



Coulomb corrections are taken into account in the simple energy gain/loss method, using a slightly higher incident 
and scattered electron energies at the vertex than measured in the lab. The shifts are the same in both cases, and 
are given by [8| : 



where Rq (in units of GeV) is given by: 



V = 0.775 - a F {Z -1)/R 



R = l.lA^/^+0.86A^ 1 /^ 



(12) 



(13) 



VII. DISCUSSION 

To illustrate the main features of the fit, the response function F^ is plotted versus W for He (left) and Fe (right) 
in Fig. [2] The top curve in each plot is the sum of the there components: a quasi-elastic peak centered on M, a 
smaller but broader "MEC" peak centered near W — 1.05 GeV, and the inelastic continuum. The main two features 
that can be noticed are that the quasi-elastic and A(1232) peaks are wider in Fe than in He (due to the larger Fermi 
momentum), and the MEC contribution is relatively larger in Fe than in He. 
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FIG. 2: Illustration of the VK-dependence of the structure function F 2 {W,Q 2 ) for He (left) and Fe (right), for Q 2 = 0.6 GeV 2 . 
The upper curves are the sum of the three lower curves (quasi-elastic, "MEC", and inelastic, from left to right). 



The Q 2 dependence of F\ is illustrated for C and Al in Fig. [3] Note the very prominent quasi-elastic and A(1231) 
peaks at low Q 2 , which "disappear" rapidly at higher Q 2 . This feature is what makes at fit accurate to better than 
10% so difficult for Q 2 < 0.2 GeV 2 . 

The fit is compared with world data [H0I ^ or ^ e ' ^> ^' anc ^ Fe/Cu in Figs. HE For W > 1.2 GeV, the agreement is 
generally within 5%, and better than 3% on average. In the quasi-elastic peak region, some larger oscillations around 
unity can be seen due to the difficulty in matching the precise shape of the quasi-elastic peak with actual data, but 
on average the agreement is within 5%. The biggest discrepancies are seen at low Q 2 , where the limitations of the 
plane- wave impulse approximation for quasi-elastic scattering become most apparent. Figure [8] shows the frequency 
distribution of the deviations (in percent) between data and fit for all the data points used. 





W (GeV) 

FIG. 3: Illustration of the IT-dependence of the structure function Fi(W,Q ) (per nucleon) for C and Al for a wide range of 
Q 2 values (units are GeV 2 ). 

VIII. SOURCE CODE 

The stand alone FORTRAN source code (named F1F209I) for this fit is available @. 
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FIG. 4: Ratios between fit and world data for He in six Q bins as a function of W. 
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FIG. 5: Ratios between fit and world data for C in six Q bins as a function of W . 
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FIG. 6: Ratios between fit and world data for Al in six Q bins as a function of W. 
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FIG. 7: Ratios between fit and world data for Fe and Cu in six Q bins as a function of W. 




FIG. 8: Frequency distribution of deviations between data and fit. 



